Single-experiment-detectable multipartite entanglement witness for ensemble 

quantum computing 
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In this paper we provide an operational method to detect multipartite entanglement in ensemble- 
based quantum computing. This method is based on the concept of the entanglement witness. 
We decompose the entanglement witness for each class of multipartite entanglement into nonlocal 
operations in addition to local measurements. Individual single- qubit measurements are performed 
simultaneously; hence complete detection of entanglement is performed in a single- run experiment. 
In this sense, our scheme is superior to the generally used entanglement witnesses that require a 
number of experiments and preparation of copies of quantum state for detection of entanglement. 
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I. INTRODUCTION 

In view of the rapid development of the experimental 
realization of quantum information processing (QIP) by 
using different physical systems, it is an urgent obligation 
to find a proper strategy for detection of entanglement. 
One of the most commonly used systems in the study of 
physical realizations of QIP is an ensemble system and 
the best-known technology therein is nuclear magnetic 
resonance (NMR) [J. 

The ensemble system of NMR has been employed for 
implementation of even relatively complicated quantum 
algorithms 0, d, 0|- However, NMR has been facing 
difficulty with regard to the existence of entanglement 
[H, [f| . The complication actually arises because of some- 
what confusing ensemble behavior of the NMR system. 
The confusion manifests itself when macroscopic quan- 
tities, such as results of NMR measurement involving 
an average over a large number of molecules, are used 
to detect microscopic properties, such as entanglement 
between pairs of nuclear spins inside each molecule. In 
fact, it has been shown that, for some particular case of 
NMR implementation of quantum nonlocal algorithms, 
the apparent nonlocal behavior of the highly mixed states 
in NMR is due to a large number of molecules involved 
in the ensemble system Q. Hence, highly mixed states 
of NMR are separable and cannot be used for immacu- 
late implementation of quantum nonlocal algorithms for 
which entanglement is believed to play an essential pre- 
requisite role d, @] . 

On the other hand, ensemble quantum computing 
(QC) pertains some inevitable advantages. It is particu- 
larly workable since spin manipulation is performed easily 
by applying corresponding pulses. In addition, ensemble 
QC, such as NMR, is supported by a long term research 
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in the area of spectroscopy. Thus, it would be unfair and 
also inefficient to totally ignore ensemble QC. We have 
been involved in realization of QIP by means of electron 
nuclear double resonance spectroscopy (ENDOR) [iol |. 
Although this system should be taken as an ensemble 
system at the moment, it has been profoundly evaluated 
for more elaborated nonlocal QIP through experimental 
studies [HI, [13]. However, after applying entangling op- 
erations in an experiment, it might still be premature to 
claim that the state is properly entangled enough for im- 
plementation of nonlocal QIP. At least some qualitative 
detection of general multipartite entanglement for a par- 
ticular system of interest should be examined in advance. 

Detection of a general multipartite entanglement is 
one of the most challenging problems for an experimen- 
tal study in QIP. There are several approaches in this 
context. One may first employ a full state tomography 
with which the complete density matrix would be ob- 
tained. Then, direct application of an existing entangle- 
ment measure would be evaluated on the quantum state 
in order to extract information about the entanglement 
of the state. This is, however, a bit too general to be 
efficient. The density matrix of the quantum state in- 
cludes far more information than necessary for an en- 
tanglement estimation only. Furthermore, it is difficult 
to find a sufficient condition of entanglement for a given 
state in general if the dimension of the state is neither 
2x2 nor 2x3 for which the Peres-Horodecki criterion 
[l3l [3] would be applicable. Detection of entanglement 
of a quantum state through violation of the Bell inequal- 
ities also should not be taken to be perfect since there 
are entangled states that do not violate any known Bell 
inequalities [HI]. 

Still one may try the existing approaches that work for 
a state that is totally unknown in advance [HI, [l?], EH • 
However, this may not be the most proper choice for our 
system of interest since as long as we are working in an 
experiment, we have definitely some knowledge about the 
state, i.e., through the prepared initial state or applied 
operations. If this is the case, it would be more appro- 
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priate to define an entanglement detector measure that 
is easily workable with less experimental effort by taking 
advantage of the available information on the state. We 
have studied the concept of entanglement witness (EW) 
[l^ . |20| with the motivation to introduce an entangle- 
ment detection applicable for a particular system of in- 
terest. The EW is an observable which has non-negative 
expectation values for all separable states. Therefore, de- 
tection of a negative value indicates the entanglement of 
the state. 

In addition to the very fast development in the th eory 
of EWs for different classes of states [H HI [H, [H , 
in experiments also detection of entanglement by the 
use of witness operators has attracted special atten- 
tion [H, H| H3, IH, IH [H- For physical systems in 
thermodynamical limits, witness operators are developed 
[3li HI, [H3, H [H, [3i|. In addition, EWs are generated 
for detecting entanglement of mixed states that are close 
to a given pure state [3?], HI] . EWs are used for charac- 
terizing different classes of a multipartite quantum state 

After determination of the EW for a particular ex- 
pected state in an experiment, the important task is to 
decompose it into local operators that are easily mea- 
surable in the given physical system. This approach for 
detection of entanglement is operationally possible and, 
for most cases, is a simple method [1^, H3, El . In Refs. 
[42l |43| detection of multipartite entanglement with few 
local measurements is studied. Throughout the study on 
this issue, it has been more appreciated to operationally 
simplify the entanglement detection process by modify- 
ing the required observables for the particular working 
system and/or decreasing the number of local measure- 
ments. 

Generally speaking, previously introduced methods for 
detection of entanglement by the use of EWs require sev- 
eral projective measurements on the copies of a state in 
order to extract the outcome. However, it is somehow 
operationally difficult to prepare several copies of a quan- 
tum state to be measured for detection of entanglement. 
Therefore, it is more advantageous if an EW works just 
with a single run experiment. Also, the detection process 
would be better to be specially modified for the particu- 
lar working system. For instance, if the physical system 
of interest is an ensemble system, then the available mea- 
surements are ensemble average measurements. 

We work with ensemble systems. Therefore, ensemble 
average measurements should be employed here, in con- 
trast with the projective measurements widely used in 
the previous works. Therefore, in this work, we propose 
a proper single-experiment-detectable (SED) EW for the 
ensemble system QC. The most significant result of our 
work is that our schemes require only a single-run ex- 
periment to detect entanglement by using nondestructive 
ensemble average measurements, which allow us to mea- 
sure several non-commuting operators simultaneously. A 
well known example of nondestructive measurement is 
the free induction decay measurement that is frequently 



used in NMR quantum computing. 

This paper is organized as follows. First in the follow- 
ing section we will briefly summarize the concept of the 
EW and will give a short review particularly for multipar- 
tite EWs. Then, in the third section we will introduce 
and prove a method with which the conventional EW 
can be transformed into a collection of separate but si- 
multaneous measurements of individual local systems of 
the ensemble QC. The analysis of this section is based 
on the assumption that the density matrix is diagonal 
after application of the disentangling operations intro- 
duced there. In the subsequent section, a different SED 
EW, which works without any assumption on the den- 
sity matrix, is introduced. As a drawback, however, we 
have to introduce an ancillary qubit. In Sec. [V] we will 
extend discussions on our scheme, study the complexity 
of the corresponding quantum circuits for SED EW, and 
give some remarks on the behavior of this scheme under 
generally existing noise in the system. 

II. EW FOR OPERATIONAL DETECTION OF 
ENTANGLEMENT 

A density matrix p is entangled if and only if there 
exists a Hermitian operator W = W\ called an entan- 
glement witness (EW), such that 

(Tr(Wp)<0, 

\ Tr(WV) > 0, Vd e S, 1 > 

where S denotes the set of separable states 0, |45| . 
Therefore, it is concluded that a state is entangled if 
some negative value is obtained in measurement of an 
EW. 

There are several methods to construct an EW [io| . 
For the case in which a density matrix has a negative 
eigenvalue when partially transposed, the construction 
of the EW is very simple. The partially transposed pro- 
jector onto the eigenvector corresponding to the negative 
eigenvalue of the partial transpose of the state is an EW 

One important and notable point about EWs appears 
when dealing with the multipartite case. EWs can be 
used to detect different kinds of multipartite entangle- 
ment by defining the space S in Eq. JT]) to be the set 
of states that does not have a special kind of entan- 
glement to be detected [39j . Let us give an example. 
While all the entanglements for two qubits are actu- 
ally equivalent to each other, for three qubits there are 
two classes of genuine pure tripartite entangled states 
[4^ |. The entangled states from different classes can- 
not be transformed into each other by local operations 
and classical communications (LOCC). One class is the 
Greenberger-Horne-Zeilingcr (GHZ) class that includes 
entangled states LOCC-equivalent to [4l| 

|GHZ) = 4=(|000) + |111)). (2) 
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The other tripartite class of entangled states is the 
Werner (W) class whose representative vector is 

|W} = -^(|100> + |010> + |001>). (3) 

This classification can be extended to more general tri- 
partite mixed states. We remind the reader that a tripar- 
tite (and generally any multipartite) state p is separable 
if it can be written as a convex combination of fully sep- 
arable states. However, a state p may not be fully sepa- 
rable but biseparable, i.e., it can be written as a convex 
combination of biseparable pure states. If a tripartite 
state p is neither fully separable nor biseparable then it 
is entangled, in either the GHZ class or the W class [37j • 
Corresponding to each class of entangled states, there are 
EWs already known. Needless to say, these entanglement 
witnesses are not sufficient to decide if the system pos- 
sesses a genuine multipartite entanglement. For instance, 
for the GHZ state, the EW is 



V^ghz = -1- | GHZ) (GHZ |, 



while for the W state, it is 



WV = ^1-|W)(W|, 



(4) 



(5) 



where 1 is the identity operator [371 ] . 

The EW that detects a genuine tripartite entanglement 
of a pure state \*S?), and of states close to is given by 



where 



w = ct - |*)(*| 



max 



(6) 



(7) 



and B denotes the set of biseparable states. Then 
Tt(Wpb) > for all biseparable states pb and 
Tr(W|*)(*|) < 0. The coefficient c in Eq. Q is de- 
termined by using the Schmidt decomposition [421 ] . 

For an ensemble system, it is desirable to find an ob- 
servable that can be performed by a small number of ex- 
periments to satisfy operational requirements. Any EW 
that is defined for the state \tp), which is the most ex- 
pected state after the applied operations and is close to 
the experimentally realized state p, should be decom- 
posed into a linear combination of individual polariza- 
tion operators, in order to make it locally measurable. 
In this paper, we introduce a new method with which an 
EW for multipartite states can be measured just with a 
single-run experiment. 



III. SINGLE-EXPERIMENT-DETECTABLE 
ENTANGLEMENT WITNESS 

For starting up an experiment, the state is supposed 
to be initialized to a simple fiducial state, such as \ip) = 



|0i...0„), for an n-qubit system. From now on we will 
drop the subscripts 1,. .. ,n to simplify notation, unless 
otherwise stated. However, in an ensemble QC, the initial 
state is prepared in the form of a pseudopure state as 
follows 



p=(l-e)l/2" + eMM, 



(8) 



where e characterizes the fraction of the state It is 
important to note here that the confusion regarding the 
concept of entanglement in an ensemble system is intrin- 
sically apart from the concept of pseudopure state as an 
initial state. A pseudopure state is used for making the 
required input state for QC. Improvement of experimen- 
tal conditions above some definite threshold is required 
to realize an entangled state experimentally. In other 
words, the pseudopure state still can be used for realiza- 
tion of a genuine entangled state if the experimentally 
required conditions are all satisfied. Here, we do not en- 
ter into the discussion of the very large required number 
of steps for making a pseudopure state as we suppose 
that, anyhow, we are given some prepared input state 
for which the status of entanglement should be studied. 

In order to produce a particular entangled state, the 
corresponding entangling operation V is applied on a 
state \ip) to yield 

Win = V\i/>). (9) 

The corresponding pseudopure state is 

Pin = (l-e)l/2 n + e|^)i nin <V|. (10) 

Note that, even though the entangling operation V is 
applied, p- m may or may not be entangled, and our task is 
to detect entanglement of p ln in order to examine whether 
p- m is applicable for some quantum nonlocal algorithm, 
for example. 

The corresponding conventional EW, Eq. ([6]), for |^)i n 

is s m s 



cl - \tp) 



= cl-V\ip)(i]j\V\ (11) 



where c = c(\ip)- m ) is determined properly so that we do 
not get a negative value for separable states [42j. The 
entanglement witness W conv detects p- ln as entangled if 
e > eiimit where 



eiimit 



Tr(W/ conv ) 



Tr(iy c0 nv) 



1 i 



(12) 



For instance, if a state pghz is the experimentally 
achieved pseudopure state for the state |GHZ) Eq. f2}, 
then eiimit is 5/7. This large value of e is in accordance 
with the results of studies on entanglement of pseudopure 
states which are indeed mixed states One may still 
work with entanglement for an ensemble system without 
being engaged in the concept of a pseudopure state pTl ]. 
However, if the state should be used for a quantum com- 
putation it is a formidable task to prepare a pseudopure 
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state that works exactly as a pure state under unitary 
operations for QC. 

Generally speaking, the expectation value of the ob- 
servable Wconv is given after several measurements on 
copies of state p ln . If Tr(pmW conv ) < then p ln is entan- 
gled. In this work, we introduce a strategy with which 
the EW can be detected in a single run- measurement. 

The most usual and convenient measurement in an en- 
semble QIP is the spin magnetization Zi of the ith spin. 
Then, we find the proper n-partite EW, which we call 
W, that satisfies the following equation 



Tr (p in W) = Tr (p in W conv ) , 



(13) 



and is detectable in a single run measurement as long as 
W is written as 



W = a + C/ f 



' a k I 



(g)n—k 



U. (14) 



The unitary operator U should be appropriately defined 
in addition to the coefficients {ai}™ =0 . It should be noted 
that Eq. ([13]) is state dependent. 

Now we will show that there exists some unitary oper- 
ator U in addition to the set of coefficients {ai}" =0 such 
that Eq. (Ti~3|) is satisfied. We also give an explicit ex- 
ample of the solution. The unitary operator U includes 
the inverse entangling operation V^ 1 and V' that is 
introduced for Eq. I|13[) to be satisfied as in Fig. [TJ 



u 



V 



out 



V 



<z> 



FIG. 1: It is possible to use a unitary transformation U and 
individual polarizations of output qubits to find the value of 
Tr(pi n Wconv). It is assumed that p ou t is diagonal. See the 
text for details. 



Suppose Pout) the state after applying the inverse en- 
tangling operation V~ x , is a diagonal matrix with a pos- 
sible classical correlation. Although this assumption is 
often satisfied, it can be dropped for a general proof if 
an ancillary qubit is added to the quantum circuit. In- 
deed, for this case the only required measurement would 
be on the spin magnetization of the ancillary qubit. This 
is discussed in the next section. 

In the present case, out task is to find some proper 
unitary transformation U and coefficients {afc} n =0 such 
that Eq. (TIB")) holds under the condition that V ^ pi n V is 
diagonal. We first prove by induction that Eq. (fl~3l) is 
satisfied for any n and later give an explicit example of 



tripartite states. By setting U — V and 



we find for the left and right-hand sides (LHS and RHS) 
of Eq. CG 



(LHS) 
(RHS) 



Tr(p out |0...0)(0. 
Ti[p out V'(b 

n 
k=l 



.0|) 



z®r- 



(15) 



1^1(16) 



Consider the case where n — 2 and let us 
denote the matrix V' by V 2 . By inspecting 
Eqs. (fTS"]) and (HHJ), we immediately notice that there 
is a solution in which the (0 . . . 0, . . . 0) element of 
V (b + YZ=i a-kl® n ~ k ®Z® I^- 1 ) V /f is -1 while all 
the other diagonal elements vanish. Note that the off- 
diagonal elements are arbitrary thanks to the assumed 
diagonal form of p ut- Typically, we can take b = —1/4, 
di = ci2 = 3/8, and 



V'o 



( o 

1 

%/3 

. J_ -i27T/3 

\ V3 




1 

vs t 

1 p -i27r/3 

1 p i27r/3 
V3 6 




1 

l 

± 

V3 



1\ 



Then, 



V' 2 (b + a x l®Z- 

I -1 






V o 



a 2 Zd 




l p i27r/3 
4 

1 „-i27r/3 



)I)V'\ 


1 -i27r/3 



i„i27r/3 





l p i27T/3 



(17) 



-i2ir/3 





(18) 



as promised. Therefore, the equality Eq. (|T3j) is satisfied 
for the case n — 2. We will take advantage of this typical 
V 2 and the corresponding parameters in the following 
proof. 

Now we use mathematical induction. Suppose that 
V n transforms 



b + a\I < 



I ® Z + . . . + a n Z <g> I ( 



I 



a Q 



to diag(— li, O2, ■ ■ • , 02") + off-diagonal terms. Then, we 
show that with a class of parameters, the unitary trans- 
formation V'n+i = Ubdn+iU P n+i{I ® V' n ) maps 

-(b + a 1 I®...®I®Z + ... + a n I ® Z ® / ® . . . <g> 7) 

+ a n +iZ® J® ...(g) I 

to diag(— li, O2, . . . ,0 2 »+i) + off-diagonal terms. The 
unitary operator U pn +i is a permutation opera- 
tion that comprises transpositions: (2,2" + 1), 
(4,2™ + 3), (2™,2™ +1 - 1). This set of trans- 

positions may be written in terms of ket vectors in 
the binary system: flOjOa . . .0„l n+ i), |li0 2 . . .0„0„+i)), 

(IdOa . . . o^a^+o, \h0 2 . . . n - 1 l n n+1 )), 

(jOila . . . l„l„+i), I li 1 2 • • • l«0 n+ i)), or equivalently 
Up n+ i = 2J \in+iin---i2h)(hin---hin+i\- (19) 

ifcG{0,l} 



5 



Thus this is nothing but a SWAP operation between 
the first qubit and the (n + l)th qubit. We also de- 
fined C/bdn+i = diag(7, H. . . . , H). This unitary trans- 
formation is easy to implement. By noting the equation 
C7 bdn+1 = (I®" <g> H) diag(ff, I,... I) we find that a sin- 
gle Hadamard gate acting on the (n + l)th qubit and 
a single C n H gate with zero-conditional-control qubits 
1, . . . , n and the target qubit n + 1 do the job. 

Now we give the details of the proof. Suppose that 
for n qubits we have the unitary operator V' n and coef- 
ficients b and {aj}™ =1 that satisfy Eq. (fT5|) so that 



An 



V' n {b + aiI®...®I®Z 
+ ... + a n Z<g>I(g>.. 



' I)V'l 



/-111 
1 



? ? 

V ? ? 



? 



? 

o ? 
? o/ 



(20) 



Here also the off-diagonal terms are not important. Then, 
for n + 1 qubits, we have the equation 



n+l 



rl 



(b + a x I( 



-a n+ iZ ( 



- a n I ( 
/(g) . 



) Z (g>I( 



A„ 
A„ 

-a„+i diag(li, . . . , 1 2 », -l2»+i, ■ 



(21) 



-1 



2 »+l, 



This should be transformed to the form: 
diag(— 1, 0, . . . , 0) + off-diagonal terms in order for 
Eq. (fT3|) to be satisfied. To this end, we apply some 
unitary transformations to B n+ \. First, we use the 
permutation f7 pn +i- Then we have 



and, for I = 2"" 1 + 1. 



,2" 



(21 - 1, 21) < — ► (21 - 2" + 1, 21), 
(21, 21 - 1) < — ► (2/, 21 - 2" + 1). 
Particularly, the off-diagonal elements of the 2 x 2 di- 
agonal blocks of B n+ i are replaced with elements that 
are 0. Now, the 2x2 diagonal blocks in B' n+ i are 
—1/2 + a n +iZ, a n +\Z, . . . , a n+ iZ , from the upper left to 
the lower right. Then we set a„+i = —1/2 and apply 
Uhdn+i — diag(J, H, . . . , H) to B' n+ i. All the diagonal 
blocks but the first one transform to X since HZH = X. 
Then we obtain 



^bd«+i-B', i+ iC/^ d „ +1 
= diag(-l,0,0,. 



, 0, 0) + (off-diagonal terms). (23) 



Therefore, the unitary transformation = 
Ubdn+iU pn+ i(I ® V' n ) and coefficients b' — 6/2, {a\ = 
ai/2}™ =1 , and a n+1 = —1/2 satisfy the equality Eq. (Tl3|) 
for all n > 2 under the imposed condition. 

The following example with n — 3 will clarify the above 
proof. For if = \(b+a l I®I®Z+a 2 I®Z®I)+a 3 Z®I®I , 
we have 



(I®V'2)K(I®V'l) 
f- 1 



2 



I e i27r / 3 
y o iV^/ 3 \e 





■i2 


i2n/3 



l„-i2ir/3 l„i27r/3 
4 e 4 e 



-i2tt/3 



+a 3 Z<g> Z® / 
-diag(-l, 0, 0, 0, -1, 0, 0, 0) + off-diagonal terms 



-« 3 diag(l, 1,1, 1,-1, -1,-1,-1) 



(24) 



Furthermore, by applying the permutation operator U p s, 
we obtain 



B'n+l — Upn+lBn+iU 



+l u pn+l 



= -diag(-l, -1,0,0,. 



. , 0, 0) + off-diagonal terms 



_idiag(l, -1, 1, -1, . 



-1).(22) 



U p3 (I®V' 2 )K(I®V'l)Ul 3 

= - diag(-l, -1, 0, 0, 0, 0, 0, 0) + off-diagonal terms 



+a 3 diag(l,-l, 1,-1, 1,-1, 1,-1). 



(25) 



It should be noted that the off-diagonal terms of Eq. 
(f2"Tj) have been also transformed under the application of 
[/p n +i. Under this transformation, the following permu- 
tations of the matrix elements, in decimal notation, take 
place. For I = 1,2, ... , 2 n ~ 1 



(21- 1,21) 
(21,21 - 1) 



(21- l,2 n + 2Z- 1), 
(2 n + 2l- 1,21- 1), 



This permutation is realized by two transpositions (2, 5) 
and (4, 7) with respect to the row labels (numbered from 
1 to 8), and the same transpositions with respect to the 
column labels. With ket and bra labels, these are trans- 
positions (|001), 1 100)), (|011),|110)), ((001|,(100|), and 
( (01 1 1 , ( 1 1 1 ) . Therefore, the off-diagonal terms, that is 
the second term in the right-hand side of Eq. 1|25|). are 
written as the matrix 
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1 „-i2ir/3 
4 e 
















1 „i27r/3 
4 e 







l„-i27r/3 
4 e 


L p -i2ir/3 















1 „i2n/3 
4 e 








1 „-i27r/3 
4 e 













l„i27r/3 
4 e 


l„i27r/3 





1 _-i27r/3 
4 e 
















1 „i2tt/3 
4 e 







1 „-i27r/3 
4 e 









(26) 



Note that all the off-diagonal elements of the 2x2 diag- 
onal blocks of this matrix have disappeared. Adding up 
the first and third terms of the right-hand side of Eq. ([23]) 
and substituting a 3 = —1/2, all the diagonal blocks are 
mapped to —Z/2, except the first diagonal block that re- 
mains as— 1/2 ( ^ ^ |. Then we use the block-diagonal 
1 \ J 

unitary transformation C4>d3 = diag(/, H, H, H) where H 
is the Hadamard operation. Then this leads to 

U hd3 U p3 (I ® V' 2 )K(I ® V'l)ul 3 Ul d3 

= diag(-l, 0, 0, 0, 0, 0, 0, 0) + (off-diagonal terms)(27) 



/7l0x0l+(l-p)llxll 



e — <z> 



Thus we found that the unitary transformation V 3 = 
U h d 3 U p3 {I®V' 2 ) satisfies V'^{\b+\aiI»I®Z+\a 2 I® 
Z® I + a 3 Z ® I ® I)V'l = diag(-l, 0, 0, 0, 0, 0, 0, 0) + 
(off-diagonal terms). This shows that the equality Eq. 
(fl3l) holds for n = 3. 



IV. 



ANCILLARY QUBIT FOR GENERALLY 
WITNESSING ENTANGLEMENT 



In this section we will show how to detect multipar- 
tite entanglement for a general case, without imposing 
the condition that the state density matrix be diago- 
nal after the disentangling operation [48|. For this idea, 
we use a single uninitialized ancillary qubit. Consider 
a single ancillary qubit initially in a thermal equilib- 
rium state: p\^ = _p|0)(0| + (1 — For a general 
proof (gen), consider an entangled state generated from 
1 ... 0) (0 ... 1 by using some entangling unitary opera- 
tion V, i.e., p gcn = |Vgo„>(^gc„| = V\0...0)(Q...0\Vl 
Then the EW is W gcn = cl — p gcn - For instance for the 
bipartite case, W gcn is non-negative for all positive par- 
tial transpose states when cd ^g-en )) is set to the largest 



Schmidt coefficient of \ip 



gen/ 



]. Our interest is not 
specific to bipartite EWs but applicable to general en- 
tanglement. 

In order to test W gen in a single NMR experiment, 
we introduce the quantum circuit in Fig. [5] Recall the 
relation 



V 



FIG. 2: Quantum circuit with a polarization measurement to 
implement the EW Wg en . 



Let us write P(0 . . . 0) = Tt(V^p ta V\0 . . . 0)(0 . . . 0|). In 
the quantum circuit, the diagonal element before the po- 
larization measurement is 

Diag {C"NOT[(p|0)(0| 

(-"NOT*}. (29) 

Here Diag A picks up the diagonal elements from a ma- 
trix A. After setting P(k) = Tr(y^ p- m V\k){k\), it is 
rewritten as 



Tr( Pin W g , 



Tr(yW^|0...0)(0...0|). (28) 



C"NOT[(p|i 

+ (l-p)|l)(l|)<g ) ^P(£:)|fc)(fc|]C n NOT t 

k 

= P(0 . . . 0)[p|l)(l| + (1 - p)|0) (0|] ® |0 . . .0)(0 . . .0| 
+ b|0>(0| + (1 -p)\l)(l\] ® £P(*)|*><*|. (30) 

The diagonal part Diag p ^J t of the reduced density oper- 
ator of the ancillary qubit before the measurement is 

P(0...0)[p|l)(l| + (l-p)|0)(0|] 

+ [l-P(0...0)]b|0)(0| + (l-p)|l)(l|] 

= [P(0 . . . 0)(1 - p) +p(l - P(0 . . . 0))]|0) (0| 

+ [P(0 . . . 0)p + (1 - p)(l - P(0 . . . 0))]|1)(1|, (31) 

where use has been made of the identity J^k = 
Trpi n = 1. Thus we have 

T±p { ±Z = (1 - 2p)[2P(0 ... 0) - 1]. (32) 
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This leads to 



P(0. 



Consequently, 



1 1 

' ~ 2 ~ 2(2p- 1) 



(33) 



TrpinWgcn = C - i 



1 



-Trp^Z. (34) 



2 2(2p- 1 

Thus, the value of Trpi n Wg Cn can be found by using the 
initial polarization TrpH Z — 2p — 1 and the output po- 
larization Tr p^^Z of the ancillary qubit. 

Extending this method, it is easy to construct a quan- 
tum circuit of concatenated EWs to discriminate differ- 
ent types of multipartite entanglement from each other. 
Suppose we want to measure k types of multipartite en- 
tanglement with a set of EWs {Wi}* =1 = {c(\tpi)) — 
F i |0...0)(0...0|F i t }fL 1 , where |^) = VS|0...0). Then 
we can test all of W% using k measurements as shown in 
Fig.H 

For example, tripartite entanglement can be differently 
written in the form of the GHZ state or the W state. For 
a general tripartite state p this can be checked by apply- 
ing disentangling operations for the GHZ state (V'ghz) 
and the W state (V^ ). This is shown in Fig. [3] if sub- 
stitutions are made as k = 2 and V\ 1 = ^'ghz an d 
V'2 — y'w ■ The entanglement of the state pi n would 
be detected by measurements of the ancillary qubit. 

In order to realize the quantum circuit Fig. [3J with a 
system such as NMR, one needs to make several free- 
induction (nondecay) measurements in a single run of a 
NMR experiment. In a conventional NMR experiment, 
however, generally one free-induction decay (FID) mea- 
surement is made at a particular point where the exper- 
iment is ceased. Indeed, there have been experiments 
involving multiple FID measurements, such as the Cory- 
48 [49| pulse sequence. 

Although it is not of particular interest to NMR re- 
searchers, in principle it is possible to take an instant 
free- induction measurement. One possible way is to ap- 
ply a 7r/2 pulse and take signals of the precession for some 
microseconds duration (that is, small enough to keep co- 
herence), and then apply a — n/2 pulse. Thus, it is a re- 
alistic idea to have multiple free-induction measurements 
in a single run of an experiment. 



V. CIRCUIT COMPLEXITY AND NOISE 
BEHAVIOR 

A. Circuit complexity with diagonal p out 

The number of two-qubit quantum gates to compose 
V n is 0(n 3 ). This is clear from the quantum circuit of 
V], in Fig. Q] Let us write the number of basic gates 
to compose V'\ by G(n). Then G{n) = G{n - 1) + 
0(n 2 ), from which we obtain G(n) = 0(n 3 ). Here we 



used the fact that an in-place C n A gate (A is a 2 x 2 
unitary matrix) is composed of 0(n 2 ) two-qubit gates 
(see, e.g., Ref. [50], p. 184). If we can use a highly selective 
pulse, the gate C n H may be performed in a single step, 
although this usually takes a long time of order 0(2"). 

In addition to the circuit complexity to compose V 7 ^, 
the circuit complexity for the inverse of an entangling 
operation V should be considered as well. Usually, this 
circuit complexity is less than 0(n 3 ). For example, a 
quantum circuit to generate a GHZ-like state from some 
pure initial state can be composed of several NOT gates, 
one Hadamard gate, and n — 1 controled-not (CNOT) 
gates. Consequently the total circuit complexity is usu- 
ally Q(n 3 ). 



B. Circuit complexity in the method using an 
ancillary qubit 

In case an ancillary qubit is employed, the dominant 
circuit complexity is due to the C n NOT gate. On the 
assumption that the internal circuit of disentangling op- 
erations (here, X may be GHZ, W, etc.) have circuit 
complexities on the order of < n 2 , the total circuit com- 
plexity of the quantum circuit of Fig. [3] is 0(n 2 x k). 



C. Noise behavior 

Although the proposed methods enable a single experi- 
ment to detect entanglement without copies of states, the 
size of quantum circuits used in the measurement process 
is considerably larger than that for usual entangling oper- 
ations. Thus noise (namely, probabilistic errors) in quan- 
tum gates can skew the result of entanglement detection 
with a higher probability than conventional entanglement 
detection using multiple copies of a state, assuming that 
error during preparation of copies of a state is negligible. 

In a simple model [51| , we assume that a quantum gate 
U acting on the target block t of qubits suffers from a 
noise such that a desired unitary transformation takes 
place with success probability p s ; otherwise, a reduced 
density operator acting on t becomes a maximally mixed 
state with failure probability 1 — p s . The superoperator 
of this noise is 



£ P .,u,t(p) =Ps(i®u) P {i®u r ) 

+ (l-p s )(Tr t p)®/W/2 lcn(t) , 



(35) 



where len(t) is the number of qubits in t; jM /2 len ^' is 
the maximally mixed state of t and I in the first term is 
the identity matrix acting on the rest of the qubits. 

Although this noise model is rather simple, analyti- 
cal calculation of the output of a quantum circuit of our 
interest under this noise is quite complicated and does 
not result in a tractable equation. Instead, we present 
a result for a particular example of the quantum cir- 
cuit shown in Fig. [1] in the case of a three-qubit input 
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FIG. 3: Concatenation of the quantum circuit of Fig. [2]to test Wi's. 



v: 



n 



-¥- 



H~H 



v: 



n- 



FIG. 4: Quantum circuit to compose V' 



1®3 



starting from the thermal state [p|0)(0| + (1 — p) 1 1) (1| 
We set the success probability as p s = h for single-qubit 
gates and p s = h 2 for CNOT gates. Noise is assumed 
to exist in individual quantum gates including the entan- 
gling operation V. We investigate the output Tr(pi n W) 
in Eq. (fl~5|) under the above noise. We choose the con- 
ventional entanglement witness for the class of the GHZ 
state given in Eq. (|4]) and decomposition into SED EW 
using the method that we have introduced in Sec. 3. Nu- 
merical results of outputs of entanglement witnesses are 
plotted against p and h in Fig. \5\ As illustrated in the 




1 0.8 

FIG. 5: Output of conventional entanglement witness 
(CONV) and that of single-experiment-detectable entangle- 
ment witness (SED) without ancillary qubit as functions of 
the probability p of starting states and the success probability 
h of quantum gates. 



figure, positive values are returned in the range of p in 
which the noise- free entanglement witness (i.e., the case 



of h = 1) returns negative values. A SED EW is more 
fragile against noise than a conventional one in the sense 
that the range of h in which negative values are returned 
is small. Nevertheless, it is important that we never ob- 
tain a negative value for a separable state even under 
noise in this example. For a general mathematical 
proof for returning non-negative values for all separable 
states is not easy because it is strongly dependent on the 
structure of the quantum circuit that has been used for 
the scheme introduced in this paper. 

One way to improve the robustness is to increase the 
value of p s . Recently, a high gate fidelity was reported by 
using a classical numerical optimization in NMR quan- 
tum computing (e.g. Ref. |52[) and this technique may 
be applicable for this purpose. This technique can also 
reduce the number of pulses needed to implement a large 
quantum gate. Thus it seems a possible way to make pro- 
posed entanglement detection methods practical in the 
future. 



VI. SUMMARY 

We proposed two schemes to reconstruct an entangle- 
ment witness into nonlocal operations and local measure- 
ments so that a single experiment without copies of a 
state is sufficient. In one scheme, an ancillary qubit is not 
required but the quantum state must satisfy some condi- 
tion, while an uninitialized ancillary qubit is required in 
the other scheme where no condition is imposed on the 
quantum state. Computational complexities and noise 
behavior have been discussed. 
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